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Abstract. The Tethered Particle Motion (TPM) method has been used to observe 
and characterize a variety of protein-DNA interactions including DNA looping and tran- 
scription. TPM experiments exploit the Brownian motion of a DNA-tethered bead to 
probe biologically relevant conformational changes of the tether. In these experiments, a 
change in the extent of the bead's random motion is used as a reporter of the underlying 
macromolecular dynamics and is often deemed sufficient for TPM analysis. However, 
a complete understanding of how the motion depends on the physical properties of the 
tethered particle complex would permit more quantitative and accurate evaluation of 
TPM data. For instance, such understanding can help extract details about a looped 
complex geometry (or multiple coexisting geometries) from TPM data. To better char- 
acterize the measurement capabilities of TPM experiments involving DNA tethers, we 
have carried out a detailed calibration of TPM magnitudeas a function of DNA length 
and particle size. We also explore how experimental parameters such as acquisition 
time and exposure time affect the apparent motion of the tethered particle. We vary 
the DNA length from 200 bp to 2.6 kbp and consider particle diameters of 200, 490 and 
970 nm. We also present a systematic comparison between measured particle excursions 
and theoretical expectations, which helps clarify both the experiments and models of 
DNA conformation. 
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1. Introduction. Single molecule studies are enriching our under- 
standing of biological processes by providing a unique window on the micro- 
trajectories of individual molecules rather than their ensemble-averaged 
behavior. Many of these studies are devoted to exploring the intricacies of 
protein-DNA interactions that are central to gene regulation, DNA repli- 
cation and DNA repair. The resolution of nanometer-scale distances in- 
volved in such interactions poses a significant challenge. The emergence 
of the tethered particle motion (TPM) method offers a practical and rel- 
atively simple solution. In this method, a biopolymer is tethered between 
a stationary substrate and a micrometer-scale sphere (a "bead" ) , which is 
large enough to be imaged with conventional optical microscopy (Fig. [T]). 
The constrained Brownian motion of the bead serves as a reporter of the 
underlying macromolecular dynamics, either by observing its blurred im- 
age in a long exposure [5], or by tracking its actual trajectory in time (e.g. 
as done in [llj and the present work) . Changes in the extent of the motion 
(which we will call "excursion") reflect conformational transformations of 
the tethered molecule. Such changes may be caused by processive walking 
of RNA polymerase [El [23] , DNA looping [3 [HI Ell EH [22l [HI [2] , DNA 
hybridization [M], DNA bending [15], HoUiday junction formation [11] or 
RNA translation [^ . 

Although TPM is simple in principle, there are a variety of technical 
challenges that must be addressed for successful implementation. For exam- 
ple, sample preparation can be compromised by multiply-tethered beads, 
non-specific adsorption, transient sticking events and dissociation of the 
tether joints [TTl [191 [iTl [3l [9] . In addition, image analysis of TPM data is 
complicated by instrumental drift and the stochastic nature of the tethered 
particle's motion. Several time scales must be considered, including the to- 
tal observation time, exposure time, and the intrinsic diffusive time scale 
of the tethered particle. We will show that quantification of the spatial 
and temporal resolution of TPM measurements requires an understanding 
of how particle motion depends on tether length, particle size and other 
controllable parameters. We focus exclusively on TPM behavior in the 
absence of externally applied force (as might be applied via magnetic or 
optical tweezers). 

The aims of this article are to: (1) review how data acquisition and 
data analysis affect TPM measurements; (2) explain a practical scheme of 
data selection and quantify the fractions of typical data that are rejected 
by each of our criteria; (3) calibrate particle motion, tether length and ob- 
servation time so that subsequent TPM experiments can be quantitatively 
interpreted; and (4) discuss the physical processes that govern TPM. Cali- 
bration of the particle motion allows precise predictions of how a particular 
conformational change of the tether, such as Lac repressor induced looping 
of DNA, affects TPM. 

Some of our experimental results were outlined in [8] . Theoretical work 
leading up to the present results on TPM motion appeared in [l3l [9l [T6] . 




Fig. 1. Idea of the tethered particle motion method. Cartoons showing the 
tethered bead in the (a) absence and (b) presence of a DNA-binding protein, 
which changes the effective tether length by looping and/ or bending the DNA. 
For example, Lac repressor protein (LacR) has two binding sites, which recognize 
and bind to two specific sequences ( "operators" ) on DNA. 

For example, Segall et al. predicted effects of cfianging tlie size of tlie bead 
and tether length, which we document experimentally in the present work. 
Our results are preparatory to experimental [7] and theoretical [16j work 
on DNA looping in the lac operon system. 

2. Results and Discussion. Using differential interference contrast 
(DIG) microscopy, the projected position of several beads in a field of view 
are recorded using a CCD camera. Sub-pixel resolution position traces for 
each bead in the image is determined using a cross-correlation method [6] . 
Standard microscopy systems such as ours are limited to two spatial dimen- 
sions; tracking of three dimensions has been accomplished using evanescent 
fields or diffraction rings [J, but this involves additional calibration and 
technical challenges. Two-dimensional tracking is sufficient for the applica- 
tions we have in mind, such as DNA-looping studies. The tracked position 
of the bead is subject to slow drift, due to vibrations of the experimen- 
tal apparatus, which we removed using a first order Butterworth filter at 
0.05 Hz cutoff frequency [19]. To quantify bead excursion, wc then used 
the square root of the sum of the variances of the drift corrected particle 
position [x, y) along two orthogonal image-plane axes: 



RMS* = ^{{x-xY + {y-yf)t- (2.1) 

Here t is the time interval over which the RMS motion is measured (typ- 
ically 4 s); a; and y represent the average of x and y over time t. Eq. 12.11 
is evaluated as a sliding filter at each point along the trajectory, and per- 
mits us to capture the tether dynamics using a single scalar quantity, as 
illustrated in Fig. [7] below. The finite-sample means x, y are subtracted 
as an additional method of eliminating instrumental drift not removed by 
the Butterworth filter; in practice, this subtraction has little effect. When 
simulating the motion numerically, we will compute the same quantity as 
Eq. 12.11 in order to make an appropriate comparison. 
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2.1. Data selection criteria. Although single-particle tracking data 
can reveal detailed features of the dynamics of protein-DNA interactions, 
care must be taken to minimize experimental artifacts such as non-specific 
binding of the bead and DNA to each other and the surface, as well as 
multiple DNA attachments on the same bead. To get acceptable calibration 
data, we implemented several selection criteria called "minimum motion," 
"motion symmetry," and "uniformity." 

"Minimum motion" discards beads that cannot be differentiated from 
beads stuck to the glass substrate (or otherwise compromised in their mo- 
bility). RMS4S from a control experiment with beads but no DNA is shown 
in Fig. [2] (green line) , and is substantially smaller than that for a tethered 
bead (blue line). Data sets exhibiting average excursions, RMS4S, lower 
than 30 nm cannot be differentiated from stuck beads and are therefore 
rejected. 

"Motion symmetry" requires that a tethered particle should exhibit 
symmetric in-plane motion about its anchor point, and is calculated from 
the covariancc matrix [3l [9] : 



are the second moments of the bead's position. Here N is the number 
of video frames and Xi, x\ arc the in-planc coordinates (i.e. the position 
x, y) of the microsphere for frame k as obtained from the drift-corrected 
data. The eigenvalues (Ai, A2) of the covariance matrix indicate the squares 
of the major and minor axes corresponding to the in-plane displacement 
of the bead and arc equal for a perfectly symmetric motion. We took 
s = \/Amax/Amin < 1-1 as our acceptable threshold. Fig. [lljd,e) displays 
scatter plots for the in-plane motion of two beads to illustrate the distinc- 
tion between symmetric and asymmetric tethers. The first plot passes the 
symmetry test and would serve as a qualified tether; the second would be 
rejected. Asymmetric bead trajectories may be caused by multiple DNA 
tethers [TT] . 

"Uniformity" qualifies tethers on the basis of the consistency of their 
motion over time and eliminates beads showing non-specific binding events, 
such as binding of DNA to the bead or glass surface for short periods. To 
detect these events automatically, we refine a procedure used in Ref. [3] . We 
first divide the entire time series into 10 subsets labeled by z = 1, ... 10. In 
subset i, we calculate RMS4S over each 4 s window and then average these, 
defining Ai = (RMS4s)i. Then we define u as the standard deviation 
of {Ai,...^io}, normalized by the overall average RMS4S. Only data 




(2.2) 



where 




(2.3) 
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Fig. 2. Selection of qualified tethers. In (a-e), the dots show instantaneous 
positions after drift subtraction; the lines show RMSab- (a) Trajectory associated 
with an accepted data set (blue) and stuck bead (green), (b) Trajectory for a bead 
that passed "minimum motion" but failed the "motion symmetry" test (see (e)). 
(c) Trajectory associated with nonuniform motion caused by transient, nonspecific 
binding, seen as a downward spike between and 50 seconds, (d) xy scatter plot 
of the trajectory in (a) shows it to be symmetric, (e) Scatter plot of the motion 
in (b) shows it to be asymmetric. The DNA used in (a-e) are 1206 bp long and 
the bead size is 490 nm in diameter, (f) Distribution of bead excursions and the 
number of beads that pass successive application of the selection criteria (see text) 
f3^. Red: original data. Cyan: after application of minimal motion filter. Blue: 
after application of symmetry filter. Green: after application of uniformity filter. 
The DNA used in (f) are 901 bp long and bead size is 490 nm diameter. 
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sets with relative standard deviation u < 0.2 are accepted. For example, 
the bead shown in Fig. ^c) meets the motion and symmetry criterions; 
however, it displays a non-specific binding event at 30s. In short, our third 
criterion removes tethers with temporal inconsistency in their Brownian 
motion. 

The first two selection criteria discard tethers that are permanently de- 
fective, whereas the third eliminates time series with undesirable transient 
events. Note that if the purpose of the experiment is to identify interesting 
molecular binding events, such as those leading to DNA looping or bending, 
then the last criterion cannot be applied, because these transient events can 
appear similar to the sticking events rejected by the uniformity criterion. 
In the present work we aimed at characterizing uniform DNA tethers, so we 
enforced all three criteria. Prior to applying the selection criteria. Fig. 
displays a broad distribution in the measured RMS4S (red). Afterwards, 
~ 50% of the data are qualified and exhibit well-defined Brownian motion 
(green bars). This figure shows that the primary cause of bead rejection is 
asymmetric in-plane motion. Experimentally, beads with multiple tethers 
can be minimized by reducing the concentration of DNA. 

2.2. Acquisition Time. The drift-corrected {x,y) trajectories are 
noisy due to the stochastic Brownian motion of the particle, and are thus 
filtered using Eg. 1 2. II over a particular time window t (usually four seconds). 
Although analysis methods exist that make no use of this windowing step 
[1], nevertheless many experiments do use it, and so we investigated its effect 
on reported bead excursion. Too short a window will increase the noise, 
leading to broad peaks in the distribution of RMSt that make signals from 
differently sized tethers too difficult to distinguish. Moreover, for short t the 
bead will not adequately explore its full range of accessible configurations, 
leading to an underestimate of RMSt , as we document below. At the 
other extreme, however, too long a window will result in a loss of temporal 
resolution. 

To determine the optimum TPM window size, we recorded data for 
200 s, for several bead sizes and a wide range of tether lengths, then found 

the mean ((RMSt )) and standard deviation (stdt = y^(RMSt) - (RMSt)^) 
of the RMS-filtered trajectory for various values of window size t (see 
Fig.[3|). Here (. . .) denotes two averages: (1) over the (200 s/t) windows that 
make up each bead's time series, and (2) over nominally identical tethered 
particles with the same bead size and tether length. The DNA lengths 
varied from 199 bp to 2625 bp, and we tested beads with three different 
diameters: 200 nm, 490 nm and 970 nm. 

Fig. [3] shows the trends as we vary t, i?bcad, and tether length L. We 
first notice that for fixed i?bcad and L, each curve levels off as t — > cxj, 
giving an asymptote that is the true RMS excursion. (For short times, the 
bead has not had a chance to explore its full range of motion in any given 
window, and so each RMSt ~^ 0, and hence so does (RMSt).) To make 
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Fig. 3. Top: Average RMS excursion and bottom: standard deviation of 
RMSt as functions of window time t m Eg. \2.1\ for different bead sizes (columns) 
and lengths L of the DNA tether (colored lines). As discussed in the text, black 
dotted lines indicate "large enough" choices oft. 



the tradeoff discussed earlier, we now ask: How long must we choose the 
window time t in order to get a reliable estimate of the true excursion? 

Naively we might suppose that each video frame gives an independent 
draw from a distribution of bead positions whose RMS value we seek. In 
that case, we would expect that as soon as i/ (30 msec) becomes large, we 
would have a good estimate of the true RMS excursion. But the top row 
of Fig. [3] shows that, on the contrary, the minimum required observation 
time increases both with increasing bead radius (moving between the three 
panels) and with increasing tether length (moving between the curves on 
a given panel). Physically, the point is that successive video frames are 
not independent draws from the distribution of particle positions, because 
the particle's motion is diffusive. The diffusion time Tdifi of a particle in a 
trap increases with increasing trap radius and with increasing viscous drag 
constant for the particle, giving rise to the trends observed in the figure. 
(For a theoretical discussion see the Supplement to [l6].) 
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Similarly, the second row of graphs in Fig. [3] shows that the scatter 
between successive determinations of RMSj decreases with increasing t. 
This "sharpening" effect also explains how RMS filtering takes rather dif- 
fuse raw data (e.g. Fig. H^-c) and transforms it into a fairly well-defined 
"state" (e.g. the individual states visible in filtered traces such as Fig. [7]). 
In both rows of Fig. [3l we have drawn dotted lines to illustrate a value of 
t that is "safe" (long enough) for tether lengths up to 2600 bp. 

2.3. Calibration of motion. In order for TPM experiments to de- 
tect discrete conformational changes of biopolymers such as in DNA loop- 
ing, it is necessary to quantify how tether length affects particle motion. 
Precise calibration data also indicates the minimum detectable change in 
tether length. Sensitive measurements may also allow detection of more 
subtle changes, such as kinking of the DNA upon protein binding or mul- 
tiple loop topologies. 
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Fig. 4. RMS excursion of bead as a function of the tether length for differ- 
ent sized microspheres, for random-sequence DMAs of various lengths. Each red 
square is the average of equilibrium amplitude of RMS motion over 20 to 200 
qualified beads, which is calculated by using Eq. \2. 1\ with t — 5 s for R — 100 nm 
(bottom data set), t = 10 s for R — 245 nm (middle data set) and t — 20 s for 
R = 485 nm (top data set). Using t = is for the same data systematically un- 
derestimates the motion of larger beads (green circles). The curves are empirical 
polynomial fits to the datasets(see Table\Vj). 

To find the empirical calibration curve, we created many DNA tethers 
of varying lengths, and attached beads of three different sizes. For each 
bead size, we estimated the RMS excursion by its finite-sample estimate 
RMSt , taking t to be the lowest "safe" value as estimated in the previous 
subsection: t = 5 s, 10 s and 20 s for beads with diameters of 200 nm, 490 nm 
and 970 nm respectively, with results shown in Fig. 31 (For comparison, we 
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Table 1 

Parameters of quadratic function ax^ + bx + c obtained for fitting both the 
equilibrium motion data (red squares in Fig.^ and 4 s interval data (green circles 
in Fig. 



Time [s] 


Diameter 


a X 10"^ 


b 


c 


5 


200 nm 


-2.58 ± 0.68 


0.17 ± 0.02 


-4.5±14.8 


10 


490 nm 


-3.37 ± 0.47 


0.19 ± 0.01 


57.3± 7.2 


20 


970 nm 


-3.49 ±0.46 


0.20 ± 0.01 


109.5 ± 8.7 


4 


200 nm 


-2.60 ±0.69 


0.17 ±0.02 


-4.75±14.7 


4 


490 nm 


-3.17 ± 0.41 


0.18± 0.01 


58.05± 6.6 


4 


970 nm 


-3.31 ±0.48 


0.18 ± 0.01 


107.7 ± 8.7 



also show corresponding results with t fixed to 4 s, which deviate signifi- 
cantly from the longer observations for the larger beads.) We summarized 
all these data with polynomial fits shown in the figure and given explicitly 
in Table m 
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Fig. 5. Theoreti cal prediction of equilibrium bead excursion, following a 
method introduced in \13l Dots: Experimental values (same as red squares 
in Fig.[^. Each dot represents 20-200 different observed beads, with the given 
tether length. Each such bead was observed for about 200 s, yielding (200s)/t 
measurements of the RMS motion, which were averaged; here t = 20, 10, and 
5 s as in Fig. Each data point shown is the average of these averages; error 
bars represent the variation (standard deviation) among the beads. Curves: The- 
oretically predicted RMS motion, corrected for the blurring effect of finite shutter 
time. For each of the three bead sizes studied, two curves are shown. From top to 
bottom, each pair of curves assumes persistence length values 47 and 39 nm, re- 
spectively, a range appropriate for the solution conditions we used J21f . There are 
no fit parameters; the theoretical model uses values for bead diameter given by the 
manufacturer's specification. The bumpiness in the curves reflects the statistical 
character of the Monte Carlo algorithm that generated them. 
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2.4. Theoretical predictions. Wc also compared the experimental 
data in Fig. 2] to a mathematical simulation of the bead-tether- wall system 
(Fig. [5]). The excursion of the bead away from its attachment point on the 
microscope slide is affected by the length and stiffness of the DNA tether, 
the size of the bead, and the various interactions between the bead/wall, 
bead/tether, and wall/tether. To account for all these effects, we modified 
the Gaussian sampling Monte Carlo technique previously used in [T3l[9llH[8] 
(see [H] for details). 

Suppose first that a semiflexible polymer chain is anchored at one 
point in space, but is otherwise unconstrained. At the anchored point we 
suppose we are given a probability distribution of different possible initial 
orientations for the first chain segment. The distribution of positions and 
orientations of the other end is then a convolution of this initial distribution 
with a kernel representing a particular diffusion process (random walk) on 
the group manifold of the three-dimensional Euclidean group. 

We can numerically compute moments of this final distribution, or its 
various marginal distributions, by a Monte Carlo procedure. Idealizing the 
polymer as a chain of finite elements, each is related to its predecessor 
by a shift along the latter's 3-axis, a twist about the same axis, and some 
random bend and twist. Rather than represent the random part using Euler 
angles, a more invariant formulation is to draw a 3 x 3 generator matrix 
from a Gaussian distribution on the Lie algebra so(3), then exponentiate 
it. The Gaussian distribution is determined by a covariance matrix, which 
represents the bend and twist elasticity of the DNA, together with bend- 
twist couplings. We estimated it up to an overall rescaling factor from 
structural data on DNA, then chose the overall factor to yield a desired 
value of the persistence length of DNA. 

Turning from the idealized problem above to TPM, we sec that we 
must implement stcric constraints: One end of the DNA tether is attached 
to a wall, which the DNA may not penetrate. Moreover, the other end 
is attached to the sphere, which itself must not penetrate the wall. Nev- 
ertheless, each segment of the intervening DNA is otherwise free to bend, 
independently of its neighbors. Thus the same Monte Carlo generation 
just described continues to be valid, except that some sterically forbid- 
den chains must be discarded. Thus our computer code generated many 
simulated DNA chains and bead orientations in a Boltzmann distribution, 
applied the steric constraints |13| . and tabulated the resulting values of 
the distance from the projected bead center to the attachment point. The 
necessary calculations were coded in Mathematics and ran conveniently on 
a laptop computer. 

We chose to compare to experimental data with "safe" values of the 
window time so we simply had the code evaluate the RMS value of this 
distance. (For a procedure valid for any t, see the Supplement to [IS]-) We 
also applied a correction to this theoretical result, to account for the bead's 
motion during the rather long shutter time (see the following subsection). 
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Fig.[5]shows that an a priori calculation of the expected motion matches the 
data fairly well, with a value of persistence length consistent with others' 
experiments; there were no other fitting parameters. 
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Fig. 6. RMS bead excursion as a function of camera shutter time in millisec- 
onds. Dots: Experimental data. Each dot represents about 20 different observed 
beads, with a tether of length L — 901 bp and a 490nm diameter bead. Error bars 
were drawn using the same method as in Fig. [5l Each point has been normalized 
to the data at 1 ms to give a dimensionless quantity on the vertical axis. Curve: 
Expected correction due to finite shutter speed, calculated by the method in the text 
(E q.\2.4]l , with shutter time given on the horizontal axis (see also the Supplement 
to fWfl 

2.5. Blurring Effect. In our experiments the camera had a long 
shutter time {5t = 31 msec). During each exposure, the bead moved, creat- 
ing a blurred image whose center is not quite the same as the instantaneous 
center. This blurring effect reduces the apparent bead excursion. Suppose 
for example that the bead has a momentary excursion to a large value of x. 
Subsequently, its stretched tether will pull it inward, so that the average 
position during the video frame has a smaller value of x. We quantified this 
effect using a 901 bp DNA and a 490 nm diameter bead at 1, 5, 10, 30, 50, 
100, 200, 300 and 500msec exposures (Fig. [S]). Longer exposures indeed 
reduce the apparent RMS motion of the bead. The effect is minimal for 
exposure times smaller than 30 ms, but decreases sharply above this value. 

These effects can be considered from a theoretical perspective (see the 
Supplement to [H]). The effect of the tether on the bead may be approx- 
imated as a harmonic restoring force. If the bead starts at a distance po 
from the center, then its average position drifts inward under the influence 
of this force. Averaging that trajectory over the video frame gives a blurred 
trajectory with center at S{po)po, where the blur factor is 

5(po) = ^[l-e-^*/^=] (2.4) 
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The time constant can in principle be estimated from first principles, 
but in practice we fit it to data such as those in Fig. [51 For very small St 
we get 5 — > 1. For large St, we have S ^ 0. 

To predict the experimental data we should thus take the theoretical 
prediction and correct it by a factor of S. This correction is trivial to 
apply (comes out of the statistical averaging), because S is independent 
of pq. The curves in Fig. [6] show that a correction of the form of Eq. 12.41 
fits the data well; this correction was applied when drawing the curves in 
Fig. El 
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Fig. 7. Two typical RMSis trajectories in the presence of Lac repressor, 
showing events of loop formation and breakdown. Total length of the DNA tether 
is L = 901 bp; bead diameter is 490 nm. Operator center to center distance is 
325. 5 bp. The upper horizontal line is the expected excursion from the calibration 
curve for the full tether; the lower horizontal line is the expected excursion for a 
tether of length 901-325 bp. 



3. Applications to DNA Looping. One of the key applications of 
the tethered particle method has been its use in studying DNA looping. 
Many transcriptional regulatory motifs involve the binding of transcription 
factors that bind at more than one site simultaneously, forming a loop of 
the intervening DNA (Fig. [T|). The TPM technique has been used to ex- 
plore these problems. The calibration analysis performed here can serve 
as the basis of a more careful evaluation of DNA looping and bending by 
DNA-binding proteins, and a guide to optimize the design of subsequent 
DNA looping experiments. For example, one may ask, what is the optimal 
total DNA length and bead size needed to reliably detect a particular type 
of loop? To answer such questions, first note that Fig. [3] quantifies how 
smaller beads and shorter tethers both allow us to work with small win- 
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dow size t, while still giving the narrow peak widths necessary to resolve 
substructure in the distribution of RMSt . Fig. [3] reinforces this point and 
also quantifies how smaller beads also optimize the resolution of TPM by 
maintaining a high slope to the calibration curve over a wide range of L. 
There are limits to what can be achieved in this way. of course: Small 
beads are hard to observe, and short DNA tethers tend to collapse (due to 
surface absorption) . Our work helps the experimenter to make appropriate 
tradeoffs when designing experiments. 

Another benefit derived from the calibration curve is a better under- 
standing of the geometry of the conformational changes we have studied. 
For example, Fig. [7] clearly shows the existence of a third state, not coin- 
ciding with either of the horizontal lines naively predicted from the cali- 
bration curve [22l [TOl [7] . More detailed simulations can then shed light on 
the geometries of the two distinct looped species disclosed by TPM assays 

mm- 

4. Conclusions. The tethered particle motion method is one of the 
simplest tools for performing single-molecule experiments on DNA-protein 
complexes. In contrast to other methods involving fluorescence, TPM never 
bleaches, allowing very long observations. The central idea is to use the 
Brownian motion of a small particle tethered to a DNA molecule as a 
reporter of the underlying macromolecular dynamics of the DNA in its 
complexes with DNA-binding proteins. The point of this paper has been 
to examine the challenges that are inherent in making useful quantitative 
measurements using this method. One of the main outcomes of that effort 
has been the development of calibration curves that illustrate how tethered- 
particle excursions depend upon both bead size and tether length. 

5. Materials and Methods. 

5.1. Sample preparation. The first step in any TPM experiment 
is construction of the relevant DNA tethers with their associated reporter 
beads. Polymerase Chain Reaction (PGR) was used to amplify labeled 
DNA with two modified primers. The primers were either biotin or digoxi- 
genin labeled at the 5' ends (MWG Biotech AG, Ebersberg, Germany). The 
labels permit specific linkage of the DNA to a polystyrene microsphere or 
glass coverslip, respectively. The PGR templates were taken from lambda 
phage or modified pUG19 plasmid (sequences available upon request). The 
PGR products were purified by gel extraction (QIAquick Gel Extraction 
Kit, QIAGEN). 

Streptavidin (Bangs lab) or neutravidin (Molecular Probes) coated mi- 
crospheres of diameter 200, 490 and 970 nm served as our tethered particles. 
In contrast to the 490 and 970 nm microspheres, the 200 nm microspheres 
were fluorescent. Prior to incubation with DNA, a buffer exchange on the 
beads was performed by three cycles of centrifugation and resuspension 
in TPB buffer (20mM Tris-acetate, pH=8.0, 130mM KGl, 4mM MgGla, 
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O.lmM DTT, O.lmM EDTA, 20 /ig/ml acetylated BSA (Sigma- Aldrich), 80 
/xg/ml heparin(Sigma- Aldrich) and 3 mg/ml casein (Sigma), filtered with 
300kD MWCO polysulfone membrane (MiUipore)). This combination of 
reagents was chosen in an attempt to maximize sample yield and longevity, 
while minimizing non-specific adsorption of DNA and microspheres onto 
the coverslip. 

The second step is DNA tether assembly. Tethered particle samples 
were created inside a 20-30 fil flow cell made out of a glass slide, glass cov- 
erslip, double-sided tape and tygon tubing. The coverslip and slide were 
cleaned with 4N HCl for 24 hours and then the flow cell was constructed in 
the same manner as described by van Oijen et al. [H]. Next, the flow cham- 
ber was incubated with 20 /.ig/mg anti-digoxigenin (Sigma) in PBS buffer 
for 30 minutes, and then rinsed with 400 ^1 wash buffer (TPB buffer with 
no casein) followed by 400 /.d of TPB buffer. Microsphere-DNA complexes 
were created by incubating approximately 100 pM microspheres with 10 
pM labeled DNA in TPB buffer for at least an hour. The DNA concen- 
tration was estimated via gel band strength. The 10:1 ratio of beads to 
DNA was designed to minimize the occurrence of multiple DNA strands 
attached to a single microsphere. The tethering procedure was completed 
by introducing 50 ^\ of the microsphere-DNA complexes into the flow cell 
for four to ten minutes. Additional tethering yield could be accomplished 
by another round of incubation with fresh microsphere-DNA complexes. 
Finally, unbound microspheres were removed by flushing the chamber with 
1 mL TPB buffer. Once microspheres were introduced into the flow cell, 
tether integrity was improved by taking care to minimize flow rates within 
the sample chamber. 

5.2. Data Acquisition and Analysis. The sample is imaged on 
an inverted microscope using Differential Interference Contrast (DIG) op- 
tics and a 1.3 NA lOOx oil-objective (Olympus). The tethered particle's 
motion was captured using an Andor Ixon camera. Each pixel dimension 
corresponds to 102 nm in the sample plane. Image transfer and storage 
was either controlled through Ixon software (Andor Technology) or custom 
Matlab code (all of our Matlab acquisition and analysis code is available 
upon request). The former recorded 8-14 bits per pixel, while the latter 
captured 14 bits per pixel. However, a comparison of the capture meth- 
ods showed insignificant differences (data not shown). Care was taken to 
ensure that the image intensity exhibited broad dynamic range without 
saturation. Some data was obtained using a Matlab-based autofocus rou- 
tine that interfaced with a Prior controller. However, for acquisition times 
shorter than five minutes, the paraxial drift was small and autofocus was 
not needed. 

The first step in analyzing TPM data is to compute trajectories for 
every tethered particle. The particle's X and Y displacement as a function 
of time was extracted from the raw data using a cross-correlation tracking 
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algorithm [5]. Such raw positional data are subject to a slow drift due 
to vibrations of the experimental apparatus. A drift correction is then 
applied using high pass first-order Buttcrworth filter at cutoff frequency 
0.05 Hz [H]. 
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